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Abstract 

We work out the holographic dictionary for the three-dimensional Schrodinger spacetimes. 
The first step in our analysis involves the correct identification of the dual sources from the radial 
expansion of the bulk fields, which turns out to be surprisingly subtle. We discuss in detail the 
holographic renormalization procedure at the linearized level and holographically compute the 
two-point functions of the energy-momentum tensor and an irrelevant vector operator. Wc 
discuss the appearance of multi-trace counterterms, parametrize the scheme dependence in our 
results, identify the non-relativistic Ward identities and compare them with expectations in the 
literature. Our results lead to valuable general insights regarding holography for spacetimes 
that are not of an asymptotically AdS form. 
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1 Introduction 



In this paper we perform a detailed holographic analysis for the three-dimensional version of the 
'Schrddinger' spacetimes of [U[2] with dynamical exponent z = 2, using three-dimensional Einstein 
gravity coupled to a free massive vector field as the effective theory in the bulk. The metric 
and vector field of the background spacetime are given by: 

ds 2 = dr 2 - b 2 e Ar du 2 + 2e 2r dudv A^da? = be 2r du (1) 

with b a real parameter. For nonzero b this spacetime cannot be conformally compactified in the 
usual sense [3j|3] and therefore it is not an Asymptotically locally AdS spacetime. Correspondingly, 
the dual theory is not expected to be conformal in the UV. Instead, significant evidence has been 
obtained [U El [7] to support the fact that the dual field theory is given by a finite irrelevant 
deformation of the two-dimensional CFT which is dual to the AdS3 solution given by (H|) with 
6 = 0. Finite irrelevant deformations of field theories generally lead to non-renormalizability and 
infinitely many counterterms and therefore result in a certain amount of non-locality and also a 
significant scheme dependence in physical observables. This non-locality and scheme dependence 
will be found in the holographic analysis as well. Our analysis will ultimately lead to explicit 
expressions for (the scheme-independent part of) the two-point functions of the energy-momentum 
tensor Ty and of the vector operator Oi dual to A^. We present these two-point functions in 
section 16.51 below. The results of this paper also remove any conceptual difficulties in extending 
the dictionary to higher-point correlation functions, although technically the computations may 
become very involved. 

Our motivations for performing this analysis are the following. First, the fact that Schrddinger 
spacetime is 'close' to AdS, in the sense that we can smoothly dial b back to zero and obtain 
precisely an empty AdS spacetime, offers us a unique opportunity to study in a controlled way the 
field theory dual to a spacetime that is not of an Asymptotically locally AdS (or AlAdS) form. The 
holographic analysis of such non-AlAdS spacetimes is in many ways still in its infancy and the results 
obtained in this paper should be very useful for further investigations in this direction. In fact, 
the gradual loss of locality resulting from an irrelevant deformation is very similar what happens if 
one gradually lowers a field theory cutoff to finite values (i.e. if one works in inverse powers of A) 
and our general results should therefore also be of relevance for a possible implementation of the 
Wilsonian renormalization group flow in the bulk theory [U El [10] . 

Another motivation for our work are several discussions in the literature surrounding in par- 
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ticular the Schrodinger spacetimes with z = 2 [TT], [12l El (13 [5] or the symmetry structure 
of the dual field theory [16]. We can add the following results to this discussion. In agreement 
with general expectations, we will find below that the non-relativistic deformation preserves a sin- 
gle ('left-moving') Virasoro algebra as well as translation invariance in the other ('right-moving') 
direction, at least at the level of the two-point functions. The left-moving Virasoro algebra includes 
the left-moving ('chiral') dilatation current which confirms the results of [3 E] regarding the exact 
marginality of the deformation from a non-relativistic viewpoint. On the other hand, we do not find 
any enhancements to an infinite-dimensional symmetry algebra in the right-moving sector. This 
contradicts the expectations of [12] and, in a sense, also of [16]. More details will be presented 
below. 

Let us briefly discuss some existing results regarding holography for non-AlAdS spacetimes. 
Many of these results are valid for spacetimes that are conformally related to AdS spacetimes, 
see for example [171 [TBI [T9l [20| I2T] for geometries dual to non-conformal branes, or [221 [23] for 
results involving the cascading theories. See also [21] for a more abstract approach. Furthermore, 
several aspects of holography for Schrodinger spacetimes have already been discussed in [5] where 
in particular the 'irrelevant deformation' viewpoint of the spacetime was already worked out in 
considerable detail. Some very interesting results have also been obtained [251 ESI ETJ [281 EHl [30] 
for the 'Lifshitz' solutions of [311 132j. These spacetimes are again not AlAdS and their holographic 
analysis is complicated by very similar issues as those described here. So far, the general approach 
for the Lifshitz spacetimes has been to impose certain boundary conditions on the bulk fields 
and renormalize the on-shell action for solutions satisfying those boundary conditions. (Related 
approaches for Schrodinger spacetimes can be found in |33[ lllj.) In contrast, our approach is to 
renormalize the on-shell action for any solution satisfying the equations of motion. This allows 
us to compute finite correlation functions of all the boundary operators rather than just a subset 
of them, but it necessarily requires a perturbative approach since in the presence of an irrelevant 
deformation the asymptotics become increasingly involved at every order in perturbation theory. 
We however regard this as an unavoidable consequence of the irrelevant deformation. It would be 
very interesting to see if our methods can be extended to the Lifshitz spacetimes as well. 

The Schrodinger solutions with five bulk dimensions were shown to be solutions of type IIB 
supergravity in |34 [ I35 [ [36] with a general five-dimensional Sasaki-Einstein manifold as the compact 
part of the geometry. Furthermore, in |36j it was shown that one may consistently truncate the five- 
dimensional field content to only the five-dimensional metric, the massive vector field sourcing the 
deformation and three additional scalars. The field theory dual to these spacetimes was argued to be 
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a so-called null dipole theory, described earlier in [37] . The three-dimensional Schrodinger solutions 
of the form considered in this paper were also shown to be solutions of type IIB supergravity in 
[7J[38], where supersymmetric solutions were found as well. The general structure of these solutions 
is very similar to the five-dimensional case and one may reasonably expect that the dual field theory 
is again a certain (two-dimensional) null dipole theory. We are however not aware of a consistent 
truncation to three bulk dimensions involving only the fields considered in this paper (i.e. three- 
dimensional Einstein gravity coupled to a free massive vector field). Although this implies that our 
results cannot be directly interpreted as a valid string theory computation (yet), we do not expect 
the correlation functions obtained from for example a 'null dipole deformed D1/D5 CFT' to be 
very different from those presented below. It would of course be very interesting to investigate this 
further. 

1.1 Approach 

The main obstacle in the holographic analysis of the background (pQ) is the unfamiliar asymptotic 
behavior of the fields near the boundary r — > oo for nonzero b. As we will see in section U this un- 
familiar behavior is exacerbated when we consider linearized fluctuations on top of the background 
dl]) and leads to radial expansions of the bulk fields which are much more complicated than in the 
AdS case. On the other hand, all our solutions turn out to be analytic in b, so when we dial b back 
to zero we can smoothly return to the linear fluctuations around the AdS background. This simple 
observation is the key to understanding holography for Schrodinger spacetimes. Namely, when we 
dial b back up from zero, the perturbative expansion in b is nothing more than an expansion in an 
irrelevant deformation parameter in the dual field theory. The holographic analysis of perturba- 
tively small irrelevant deformations was worked out in detail in |39t HO] and we can use the results 
of those papers to find the dual correlation functions up to any finite order in b. In this paper we 
perform a detailed holographic analysis up to order 6 4 and use these results to conjecture the result 
of the asymptotic analysis to all orders in b. Our conjecture passes several consistency tests and 
could only be incorrect if there would be a significant structural change appearing at order 6 5 , a 
scenario we consider highly unlikely. Combining our knowledge of the linearized bulk solution to 
all orders in b with our conjecture for the all-order results of the asymptotic analysis then allows 
us to extract the two-point functions of the dual field theory for finite b. 
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1.2 Outline 

In the next section we set up the bulk problem by introducing the bulk Lagrangian, computing 
its equations of motion and discussing the linearized analysis. We also explain in more detail how 
to treat the problem perturbatively in b. We then holographically compute the relevant two-point 
functions for b = 0, so for empty AdS3, in section This gives us the basic correlators that will 
receive an infinite series of 6-dependent corrections once we switch on b. In section 0] we present 
the solution to the linearized equations of motion for finite values of b. We explain how a careful 
b — > limit leads to a proper identification of the sources in the radial expansion of the fields. The 
analysis in this section, in particular subsection 14.21 contains the most important conceptual results 
of this paper. We then continue with the holographic analysis: first to order b 1 in section [5] and 
then to order b 4 in section [6l In fact, in section [6] we also present our conjecture for the all-order 
result and compute the associated two-point functions. In the appendix to this paper we briefly 
discuss the irrelevant deformation from a field theory perspective. 



2 Setup 

We consider a three-dimensional Einstein-Proca bulk theory with the following action: 

and we set A = — 1 henceforth. We can also set 167tGtv = 1 and we need m 2 = 4 in order for the 
Schrodinger spacetime with z = 2 to be a solution. The equations of motion take the form: 

Rulu - \RG^u - Gfiv = \f^F u P + 2A ti A u - G^F^FP* + A P AP) 

2 2 8 ( 3 ) 

V% - 4A„ = 
We will always be working in a radial-axial gauge where: 

G llv dxPdx u = dr 2 + j ij dx i dx j . (4) 

In this coordinate system the extrinsic curvature of a slice of constant r is 

K ij = \% ■> ( 5 ) 

where a dot denotes a radial derivative. 

One solution to the equations of motion is empty AdS3 with a metric and vanishing vector field: 

-/ ij dx i dx j = 2e 2r dudv A p = . (6) 
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Another solution to the equations of motion is the Schrodinger spacetime: 

j ij dx i dx j = -b 2 e Ar du 2 + 2e 2r dudv A^dx^ = be 2r du . (7) 

Notice that we work in conventions where rjijdx l dx :J = 2dudv so r/ uv = 1 and lightcone indices can 
be raised and lowered without introducing factors of 2. 

When we attempt to compute the on-shell action for the above solutions we find divergences as 
we integrate all the way to the conformal boundary at r = oo. We cancel these in the usual fashion 
by adding a counterterm action. In our case we add counterterms of the form: 

,2 



S c t = / d 2 x^(-2 + AiA % ) , (8) 



which ensures that the on-shell action remains finite for both of the above solutions. Furthermore, 
the first variation of the renormalized on-shell action 

S(S + S ct ) = J d 2 x^j{-K lj + 1 ij (-l + K + ^A k A k )-A i A j )5 lij + J d 2 x^{F ir + 2A i )5A i (9) 

also vanishes on the above backgrounds with this choice of counterterms. We note that this 'ad hoc' 
definition of a counterterm action is somewhat illegal, since in principle one should find the most 
general asymptotic solution to the equations of motion (to the desired order in the fluctuations) and 
then one computes rather than postulates the counterterm action by following the usual holographic 
renormalization methods. However this is precisely what we will do below when we compute the 
counterterm action for perturbations on top of these backgrounds. Therefore, if our initial guess 
(jHJ) for the counterterm action happened to be wrong then we will automatically compute, order by 
order in the fluctuations, the corrections that bring us to the right counterterm action. Of course, 
our reason for adding ([8]) is that it provides us with a convenient starting point. 

To find the two-point functions in the dual field theory we will below perform a linearized 
analysis of small fluctuations around these backgrounds. For the metric variations we stay in the 
radial-axial gauge so 5G rr = 0, 5G r i = and we define = hij. For the vector perturbations we 
define 5A^ = and bF^ v = fuv- Henceforth jij and Ai will denote the fixed background metric 
and gauge field, respectively, and indices are raised with j lJ . In summary, then, we have: 

7ij lij + h ij -> A n + a M • ( 10 ) 

We will need to compute the on-shell action to second order in the fluctuations. Expanding ([9]) we 
find that this on-shell action takes the form: 

S = J d 2 x^(h^-h i ^+^(h k k +A k a k + ^A k h kl A l )-2A^)h lj + J d 2 x^(f ir +2a i )a i (11) 



where is shorthand for "y lk hki^^ and similarly h)t = ^ kl hki- Notice that the bulk term in the on- 
shell action for the fluctuations vanishes because both the background and the fluctuations satisfy 
the equations of motion. 

As we announced in the introduction to this paper, the main idea of the subsequent compu- 
tations will be to work perturbatively in the deformation parameter b. In this sense b is exactly 
like the coupling constant A of the <3? 3 model considered in [39]. In particular, the solution to the 
linearized equations of motion will have an expansion in b, 

a-i = aioji + a {i}i + a{2}i + • • • . N 

(12) 

hij = h{o}ij + h{i}ij + h{ 2 }ij + ■■■ 

where the curly brackets indicate the corresponding order in b. Furthermore, the background fields 
jij and Ai, which appear in the linearization of the equations of motion ([3]) and in the on-shell 
action (jlip . also have an expansion in b. These expansions however terminate at order b 2 . Again 
this is just as in the A<£ 3 model of [39] where the on-shell action also contains an explicit A<3? 3 term. 



3 Analysis for 6 = 

Our first task is to perform the holographic analysis of linearized fluctuations around empty AdS3. 
We henceforth work in lightcone momentum space with lightcone momenta q u and q v . The analysis 
in this section closely follows standard results in the holographic renormalization literature, see in 
particular |41 | 142 1 143] . We refer to |44j for an introduction to holographic renormalization. 

3.1 Linearized asymptotic solution 

The linearization of the equations of motion ([3]) around the AdS3 background ([6]) is straightforward. 
One finds that the metric and vector field fluctuations decouple and the equations are easily solved. 
For the metric we find the following solution: 

hij = e 2r h^ij + h(2)ij (13) 

where the integration constants fyo)r/ are the sources and 

(14) 

d j h h)i - dlh (2)ij = o. 



s 



Below, we will be working exclusively in momentum space and in components. In that language 
the above equations become: 

h(2)uv = -^(lv^(0)uu — ^QuQvh(p)uv + Quh(0)vv) 

1vh{2)uv - Quh(2)vv = 
Quh(2)uv - 1vh{2)uu = • 

As expected for Einstein gravity in AdS3, this solution for hij is just a linearized diffeomorphism. 

For the massive vector field we obtain the following asymptotic solution to the linearized equa- 
tion of motion: 



ft; 



a (0)i e r + 0( 2 )i + (r5( 4 )j + a (4)i )e r + ... (16) 



where a(g)j are the source terms and the first two subleading terms are given by: 

a (2)u = -7Qu(q u a(o)v - Qva(o) u ) 

4 ! (17) 



-.3 , 

4' 



a {A)u = --:Q u Qva(o)v 



and similar equations hold with u and v interchanged. As usual, the are not determined by 
the asymptotic analysis of the equations of motion alone but will be fixed by demanding regularity 
in the interior of the spacetime. 

3.2 Holographic renormalization 

Upon substitution of the asymptotic solution into the on-shell action, we find divergences as we 
send the cutoff r — > oo. These are cancelled by the counterterm action: 

Sct = ~\l d2x Vl(l^ D 7 a k + (d k a k ) 2 + ri-^U^ak + ^d k a k a^d m a m )^ (18) 
with D~ = e~ 2r D = e~ 2r 5 l i didj the Laplacian associated to the induced metric 7^. This countert- 
erm action can be found with the usual means, in particular we obtained it using the Hamilton- 
Jacobi method of [45]. We will not repeat its derivation here. Notice that the counterterm action 
is covariant up to the explicit appearance of the coordinate r, which signals the presence of a 
conformal anomaly. 

3.3 One-point functions 

After adding the counterterms, the first variation of the renormalized on-shell action becomes finite 
and takes the form: 

S(S + 5ct,o) = I d 2 x [(4aJ 4) - V (0) + in^5 fc a (0)fc )oo (0) , + + l -R[h {Q) }^)5h [0)lj \ (19) 
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with i?[/i( )] = d l d 3 /i(o)ij — ^ 4 ^*^(o)j the linearized curvature associated to fyo)ij- The one-point 
functions are then given by: 

{&) = --= L=-*_(5 + 5 ct , ) = -4a} 4) + ^n'ajo) - ^ ! 3 fc a (0)fc 

y/~9(o) oh {Q)ij W 2 

These are one-point functions in the presence of sources, so differentiating them once more with re- 
spect to the sources gives the two-point functions, see below. The contact terms in (O l ) are scheme- 
dependent, although conformal invariance constrains them to a linear combination of d l \Dd 3 d(o)j 
and \3 al Q y We will elaborate on the scheme dependence of our results in the following sections. 

3.4 Ward identities 

Equation (|20p expresses the one-point function of the energy-momentum tensor in terms of the 
normalizable modes h^)ij- In the asymptotic analysis of the equation of motion we obtained 
the constraints (|15p for the coefficients h^)ij- We may now translate these equations into Ward 
identities involving the energy-momentum tensor. For example, the first equation in (j 1 5 j) leads 
directly to the trace anomaly: 

(V)=A = 2^R[h {0) ] = ^R[h {0) ]. (21) 

Reinstating the factor of (167tGat) -1 , we find c = 3/(2Gn), as expected [Ml EH EE] • Furthermore, 
one easily finds that the last two equations in (fT5j) become the usual conservation equations: 



d^Tij) = (22) 
where we also used the first equation in (|15l) . 

3.5 Callan-Symanzik equation 



In [30] the holographic derivation of the Callan-Symanzik equation for scale transformation of the 
on-shell action was explained in complete detail. Although we will need such a detailed analysis 
at higher orders in b, at this order the analysis is familiar and essentially dates back to [471 EH] • 
Following the known prcoedures, then, we find that the Callan-Symanzik equation takes the form: 

I d 2 x(2(g {0)ij + hrouj)— ^2a {0)iT ^—)S rcn = -B . (23) 

J x on (0)ij oa {0)i / 
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Here the integrated conformal anomaly is denoted B. It originates from the manifest radial depen- 
dence of the counterterm action, more precisely 



where the partial derivative with respect to r is taken with all the bulk fields fixed. On the left-hand 
side of (|23h all the sources scale with their canonical dimension, so there are no beta functions at 
this order. 

Notice that the integral of the conformal anomaly A appearing in (12 ip should in principle be 
equal to B as given in (I24p . However in perturbation theory the two are manifestly different. This 
is firstly because the terms appearing in A are topological invariants (and we will also find total 
derivatives at higher orders below) and therefore integrate to zero in B, and secondly because the 
terms appearing in B are all second order in the sources and therefore do not show up in A, which 
in the linearized analysis is only determined to first order. In perturbation theory, then, the Callan- 
Symanzik equation and the trace Ward identity are both valuable and complementary sources of 
information. 

3.6 Two-point functions 

The above analysis relied only on the asymptotic analysis of the equation of motion, which is always 
sufficient to obtain the Ward identitites and the Callan-Symanzik equation. On the other hand, to 
compute two-point functions we need to demand regularity inside the bulk as well. For the massive 
vector field the solution to the linearized equations of motion is given in terms of Bessel functions. 
More precisely, we find the regular solution to the equation of motion to be: 




(24) 



a u = a 



(o) u QuqvKo{y / 2q u q, : 



'V 



(4e 2r + 

) + 9u°(0)«)- K 'i(V^9u* 



•v 



e- r ) (25) 



and similarly with u and v interchanged. We can expand the Bessel functions, 



Ko(z)= (-7 + ^(1 -7)z 2 + ...)+log(z/2)(l+ l -z 2 + ±z 4 + ..: 
Ki{z) = + i(2 7 - l)z + . . . ) + | log(*/2) (l + ±z 2 + -^z 4 + ... 



) 



) 



(26) 



and compare with (|16p to find an exact match for the leading terms and to also find that 



1 

a (4)« = o 



qIqv(MQuQv) ~ hi(2) - 1 + 2j)a {0)v + — a {0)u q 2 u ql . 



(27) 
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We may now substitute this expression the first equation in (|20j) and take another functional 
derivative to obtain the two-point functions. In conventions where 

Oi(x) = J ^OMY^ (Oi ■■■) = . •) , (28) 

with q ■ x = q u u + q v v and the factor of i being a consequence of the Lorentzian signature (see 
appendix B of [49]), we find that, up to contact terms: 

(O u (u,v)O u (0)) = - l - [ -f^- q lq v \n(q u q v - ie)e^ x = (29) 



2 J (2tt) 2 ^ ' v ™™ ' 27ru 4 « 2 

and of course similarly with u and v interchanged. Notice how we inserted an ie prescription in 
this correlation function. This particular ie prescription corresponds to a time-ordered two-point 
function. We will be computing exclusively time-ordered two-point functions in this paper. The ie 
prescription can in principle be found from a careful bulk analysis, see [501 [5l] for details, the end 
result of which will be precisely the insertion as we have written it0 To enhance readability we 
have omitted the corresponding ie prescription from the position-space expression, for time-ordered 
two-point functions it is easily reinstated by replacing t — > t — iet. 

For the energy-momentum tensor the normalizable coefficients h^)ij are in fact completely 
fixed by the Ward identities, which is a reflection of the fact that three-dimensional gravity has no 
propagating degrees of freedom. We find that 

Qu 1 0^ 2 

h(2) uu = — h(2)uv = t(— h(o) vv + q u qvh( ) uu - 2q u h^ uv ) (30) 
Qv 4 q v 

and therefore, using that in Lorentzian signature insertions of the energy-momentum tensor are 
defined via 

(!*...) = -4m -A- (...), (31) 



8h i3 



we find that up to contact terms, 



(T uu (u, v)T uu (0)) = -An 2 i I .f* g " . e*** = ^ - (32) 



(27r) 2 q u q v - ie u A 

and again similar results hold with u and v interchanged. The coefficient of the two-point function 
is c/2 = 3/(4Gyv) if we reinstate the factor IQttGm again, in agreement with expectations. Notice 
that we rewrote g" 1 as q u /{QuQv — Such a prescription should always be understood when we 
write inverse powers of momenta since it is the correct regularization of the singularity at q v = 
for time-ordered correlation functions. 



1 Another discussion of the real-time subtleties in Schrodinger spacetimes can be found in [52] 
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4 Asymptotic behavior for finite b 

As we emphasized above, our approach will be to treat the parameter b as infinitesimally small and 
to work perturbatively in b. In this section we however first present the solution to the linearized 
equation of motion for finite b. Apart from certain technical complications, already discussed in 
[5], the main conceptual difficulty is the proper redefinition of the integration constants to what 
will eventually be the source and vev terms in the dual field theory. We shall demonstrate that 
this redefinition can be performed by taking a careful 6 — > limit. 



4.1 Asymptotic solution 

For finite b the linearized equations of motion are rather complicated. Exact expressions for these 
equations can be found in [5], where the asymptotic solution was also presented. Here we will 
review and slightly rewrite this asymptotic solution. 

We first consider the radial component a r of the vector field. Regarding the dynamical part of 
a r , we find that we can combine the equations of motion to write a single fourth-order equation 
involving only a r of the form: 

A 2 2 i a 2r It2 2 //\ 

4q u q v a r + 4e q u q v [b q v a r - a r ) 

(33) 

+ e 4r (-16a' r - 4< + b 2 q 2 v ((8 + b 2 q 2 ) a r - 2 (2a' r + <)) + 4a' r " + <") = 

where a prime denotes a derivative with respect to r. This equation can be solved in terms of 
Bessel functions. In addition to the solution to this equation there is also a non-dynamical part of 
a r which comes from a linearized diffeomorphism. The most general solution for a r is then of the 
form: 

ib 

a r = —X{2)w + ciXi(e" 2r q u q v ) + C-il-i(e~ 2r q u q v ) + c 9 X 9 (e" 2r q u q v ) + c_ 9 X_ 9 (e" 2r q u q v ) (34) 
qv 

where X(2)ot> c ±i an d the c±g are five integration constants, which a priori can be arbitrary functions 
of q u and q v but not of r. (The notation X(2)w will be useful below.) The functions I±i(x) and 
I±g(x) take the form: 



I±i(a;) = V2x( Sl ± 1)I± s1 {V2x)tI-i± s1 W2x) 

(35) 

Zfc 9 (x) = V2x(s 9 =F 3)I ±S9 (V2x) =F I-i± S9 (V2x) 



with 



si = y/1 + b 2 q% s 9 = y/9 + Wq% (36) 



13 



and I u (z) the modified Bessel function of the first kind. Recalling the series expansion of the Bessel 
function, 

we see that the four solutions have asymptotic behavior: 

l±i(e- 2r q u q v )~e-( 1±s ^ r + ... l± 9 (e~ 2r q u q v ) ~ e^ 1 ^ + . . . (38) 

where the dots represent infinite series which trail the leading term by powers of exp(— 2r). No- 
tice the rather unusual behavior where the power of exp(r) becomes dependent on the lightcone 
momentum q v . This phenomenon is typical for the Schrodinger spacetimes and shows that the 
anomalous dimensions in the dual theory become dependent on the momenta. 

As usual, regularity inside the bulk will eventually require us to pick specific linear combina- 
tions of these four Bessel functions. More precisely, the regularity condition will lead to two linear 
relations between the four coefficients c±\ and c±g. It is however important that we only impose 
such a regularity condition after finishing the asymptotic analysis. This is because the holographic 
renormalization procedure is equivalent to finding the correct variational principle for the gravi- 
tational dynamics |53| . As in ordinary classical mechanics, obtaining such a variational principle 
requires one to treat the conjugate momenta as unrelated to the position variables. On a practical 
level, regarding the source and vev terms as independent becomes essential for obtaining correct 
results in the presence of multi-trace operator mixing }40j . For these reasons we will not yet impose 
regularity in the interior. 

Using the above solution for a r one can easily reconstruct the solution for the other component 
fields using the remaining equations of motion. For the metric components we find that: 

i _ 2r i 

vv 

1. 



,2„4r,, r.2„2r_,. , „2r , 

a 



n uv = — -o e X(o)w -be rx(2)w + e " X(o)uv + X(2)uv 



h uu = b e r — (-q u qvX{o)w + ^X(2)w) - 2b e r r— X(2)w + e r X(o)uu + X(2)uu 

qv Qv 



+ ciUi(e r q u q v ) + C-{H-i{e r q u q v ) 

where we introduced five new integration constants X{o)uu, X(o)uv, X(o)w, X(2)uu and X(2)uv, whilst 
X(2)w an d c±i were already introduced above. The highlighting will be explained later but we warn 
the reader that the colored coefficients will not be equal to the sources of the boundary theory. 
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(40) 



The functions H±i are slight variations of the combination of Bessel functions T±i, namely 

U±i{x) = T8iV2bqlq v x- 5 / 2 (6 + 2b 2 q 2 v T 6si + x)J ±ai _ 4 (V2s) 

+ 32ibq 2 u q v x-' 3 (b 2 q 2 v (* a T 6) + 12 + (« x =F 2) x) / ±sl _ 3 (\/2^) . 

Furthermore, the coefficients X(2)ij are related to the X(o)ij via: 

X(2)uv = j(QlX(0)uu - 2quqvX(0)uv + QlX(0)vv) 

QvX(2)uv - QuX(2)vv = 
QuX(2)uv - QvX(2)uu = 

which happens to be identical to (|15p in the undeformed theory. 

For later reference we will give the asymptotic expansion of h uu . Using (|37p we find that: 

h uu = b e —{-q u qvX(o)vv + 4 X(2)«J - 26 e r— X(2)w + e X(o)uu 

, K 3- S1 ) (_ ^i-l^fe^)^)/ 2 Z2§~l^(gl - l)(g^) (3+Sl)/2 -2, \ 

Cl { ^i(3 + 6 2 g2-3 Sl )r( Sl -2) &^r(i + Sl ) e 

, r(3+ «) / ^i+s^b^)^)/ 2 i2l+^H-^i-i)(gng,) (3 - sl)/2 _ 2r 
c_1 ^ ^i(3 + 6 2 (?2 + 3 Sl )r(- Sl -2) + ^r(i-si) 6 

i23 + 3«( ?ttfe )(5-i)/2 ^_ 4r+ \ 



(42) 



6(?3 Sl (-3 + 6 2 g2)r(- Sl -2) 

For a u and a v we find similar expressions as for a r which can again be written in terms of four 
Bessel functions. We will present only their asymptotic expansions here. These are given by: 

,(i-«.K/„., f -2Q»l-l) , d(~2 + »l) .-4, 



, (1+ Sl )r , A + b2 Qv + gl) -4r , 

+ e ^"i 1 " 86 2 (6 2 g 2 -3) 6 + - 

n- 89 )r,/. A , _QuQ^-2r , Qlqy(3b 2 q 2 + 5 + 5s 9 ) 4r 

+ ^ (4)U V 1 + ^9 + 8(ll + 6 2 g 2 + 3 S9 )(3 + 6 2 ( ? 2 + S9 ) e + " 

_ {1+sg)r / 2(s 9 -3) 2g M (18-6 S9 + 6 2 gg(2- g9 )) _ 2r g g(16 + 3b 2 g 2 - 8s 9 ) „ 4r 
ll q v (8 + b^)s 9 6 + 4(8 + 6 2 (? 2 )(-2 + S9 ) e 



+ e 2 ^ 



lu 26 \ 1 

— X(0)vv - -^X(2)vv I + 2 b *(2)t 



2<?„ 

(43) 
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and 



■<i+'i>',/,.., I 1 j- *^ Ha ".-" 



+ e -(i+«»)r^, f g'( s 9~ 2 ) | gugg(-2 + a 9 ) 2r \ 

V(4)« ^ 2(3 + 6 2 g 2 + sg) + 4(1 + Sg)(3 + b 2 q 2 + sg) 

1 l 2r 1 z. 

+ X(0)w + 2 b X(2)™ • 

Here we introduced new integration constants f/>(o)j an d ^(4)1 which are a related to the c±\ and 
c-1-9 as: 

23( 1 +l)« Zu (g tt(ftl )3( 1 -'0 

^ fo^o 

25( 1 +^)^(g u g j; )5(l-^) 

a x ( f , (45) 
%)u iq v (-2 + s 9 )T(l + s 9 ) C9 

^ = f(l+li) Cl • 

We will henceforth use the integration constants VYo)* an d ipu)i rather than c±i and c±g. 

4.2 Interpretation 

The solution presented in the previous subsection appears rather unstructured. In contrast to the 
6 = analysis of section [31 where every component of the various fields had a simple expansion 
in powers of exp(— 2r), here we find several expansions that appear simultaneously in various 
components of the fields. Let us now discuss how to perform the holographic analysis in this case. 

Our first task is to identify the terms in the radial expansion that correspond to the sources of 
the bulk theory. This would be a trivial exercise in the usual AdS/CFT setup, where the sources 
are usually just the leading terms in the radial expansion of the fields, at least up to well-known 
subtleties for operators with A < d/2 discussed in [53] and up to multi-trace deformations discussed 
in for example |55[ [56"1 [571 [58l 1591 [60l [6T1 HQ] . In our case inspection of the solution naturally leads 
to the idea that the V>(o)i an d X(a)ij could be sources in a dual field theory. These are the coefficients 
marked in red in the above formulas. Indeed, the ^(0)1 are leading terms in an infinite expansion 
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in powers of exp(— 2r) and upon sending b — > the power of exp(r) accompanying the coefficients 
^(0)11 an d ip(p) v becomes precisely exp(2r) in a u and a v , respectively, since in this limit si — > 1 and 
sg — > 3. The 0(o)i therefore 'land' precisely at the order of the source terms in the 6 = expansion 
of the previous section. Similarly, the X(o)ij appear precisely at the same order, exp(2r), as the 
sources in the b = case. They might therefore correspond to the energy-momentum tensor in 
the dual theory. All the other terms in the radial expansion of the fields would then merely be 
'induced' terms which are sometimes more leading than the source terms [39J . 

There is however a problem with the reasoning of the previous paragraph, because upon further 
inspection the b — > limit is not as straightforward as our previous discussion seems to suggest. 
For example, the first subleading term in the expansion multiplying ip^ v in the expansion (|43j) 
for a u also becomes of order exp(2r) as b — > and will therefore mix with to become the 

actual source term a^ u . A similar mixing can be seen to occur for many other coefficients in the 
radial expansion of the fields. The terms highlighted in red in the above formulas therefore do not 
precisely reduce to the source terms as b — > 0. 

This observation corroborates nicely with the results of [39} 140] , As we mentioned above we 
view the appearance of b as an irrelevant perturbation. In [39] it was demonstrated that to any 
finite order in the irrelevant deformation parameter, the source is always equal to the term at order 
exp(— A_r) in the radial expansion of the fields. Here A_ is the appropriate scaling dimension for 
the dual source without the deformation, which in our case happens to be equal to —2 for both the 
metric and the vector field. This simple prescription leads one directly to the correct identification 
of the sources, as we will now proceed to demonstrate. 

4.3 Identification of the sources 

The result of [39] is valid perturbatively in the irrelevant deformation parameter, i.e. in b. This 
in particular implies that the 6-dependent powers of exp(— r), so terms like exp(±sgr), should be 
viewed as infinite series expansions in b before the analysis of [39] can be applied. In order to 
identify the proper sources we should therefore write such terms as: 

exp(±s 9 r) = exp ( ± y^Tb^ r) -> exp(±3r) ± Qfc 2 ^ - ^b 4 q^ r + ^b A q A v r 2 + . . )j 

(46) 
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so as a term of order exp(±3r) plus an infinite series of logarithmic terms^ As explained in [40] . 
the presence of these logarithmic terms does not modify the identification of the source as the 
term multiplying simply the exp(— A_r) coefficient. For the identification of the sources it is then 
sufficient to replace exp(±sgr) with exp(±3r). Of course the same analysis can be applied to the 
other 6-dependent powers, so: 



exp(± Sl r) = exp ( ± y/l + Wq% r) -)■ exp(±r) ± Qft 2 ^ - I& 4 ^ r + \b\y + . . . V 



(47) 



Taking into account such an interpretation of the fractional powers of exp(r) allows us to directly 
read off the correct source terms in the linearized solutions. For example, from f)43|) and (|44p we 
directly find that the terms of order exp(2r) in the radial expansion of the are: 

2g u (18-6s9 + 6 2 ^(2-s 9 )) . bq u 2b 

O(0)« = ¥>(0)u (o + h 2 a 2\ qn no)v + IT- X(0)vv ~ -$X(2)vv 

q v {& -+- q v )Sg zq v q v 

a (o)v = ^(o)« + 2^ X (°) VV ■ 

It will be useful to also define the corresponding 'normalizable' terms, which by definition sit at 
order exp(— 2r) in the radial expansion of the fields: 

" 2 " 2 ^ (26 2 ? 2_ 5S9)(1 + Sg) 



a ^ u = ^ + 8(1 - bW, + S1 )^ U + 12(5+^)78 +7a(- 3 + s 9 ) 



2 2 

a ^ = ^ + 88 + 8bW v ~ 24s 9 ^ ~ ■ 



(49) 



For the metric we find that: 



h(0)uu 


= X(o> 




= X(o) 




= X(o) 



8b 2q u 



(50) 



are the correct sources in the dual theory, whereas the 'normalizable' terms are given by the terms 
at order one: 

Ql i 4g n ( Sl - 1) 
26(2 bfi—^* 

(51) 



h(2)uu 


= X(2)uu + 


h(2)uv 


= X{2)uv 


h(2)vv 


— X(2)vv ■ 



2 We call terms involving powers of our radial coordinate r 'logarithmic' because they are powers of logarithms 
of the conventional Fefferman- Graham radial coordinate p = exp(— r) and also because they correspond to beta 
functions and conformal anomalies in the dual field theory. 
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The above transformation between the correct field theory sources, so fyo)y an d and 
the natural integration constants, so X(o)ij an d V'(o)j) is quite remarkable. Indeed, the modes 
parametrized by the Xij represent a linearized diffeomorphism whereas the modes parametrized by 
the ipi represent a propagating fluctuation. They are thus qualitatively very different and one may 
naturally try to associate them with different operators in the dual theory. We however now observe 
that the correct field theory sources fyo)ij an d a(o)i are not in one-to-one correspondence with the 
two sets of solutions parametrized by the Xij an d the ipi. Rather, each of the field theory sources 
switches on a combination of these two solutions. This is in fact reminsicent of the analysis of TMG 
at the chiral point presented in [49], where a similar mixing between the propagating solution (i.e. 
the Bessel function) and the linearized diffeomorphism occurred. 

We would like to emphasize that the redefinitions (|48j) and fj50f) are very non-local because of the 
factor of q~ l and q~ 2 appearing in the transformation formulae. Notice that, despite appearances, 
we are using an ie prescription such that these factors are in fact non-local in the u direction, see 
for example the two-point function ([32]). Any analysis where the X(o)ij an d ^(o)i are treated as the 
field theory sources will therefore hide a non-local convolution of the correct field theory correlation 
functions and may easily lead to erroneous conclusions If| 



5 Analysis to first order in b 

Altough we have just presented the solution to the linearized equation of motion for finite values 

of b, we cannot directly use this solution to compute correlation functions since we are unable to 

perform the holographic renormalization procedure for finite b. We announced in the introduction 

that our strategy will be to work perturbatively in b instead. In this section we will make the first 

step by performing the holographic analysis to first order in b. The structure and methods of this 

section are exactly the same as that of [39\ 00] so the presentation will be rather brief. We will 

however discuss the field theory results in a little more detail. 

' 5 Although the field redefinition from the ftmjy and <i(o)i to the X(o)ij an d ip(p)i ls non-local, there do exist many 
other field redefinitions which are local and therefore in principle allowed. At large N, the full space of such local 
field redefinitions also includes local double-trace deformations and Legendre transforms akin to [54]. As explained 
further in the remaining sections, rather than exhaustively exploring all possible redefinitions we have instead opted 
to fix this freedom in a rather natural way. 
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5.1 Linearized asymptotic solution 

From the results of the previous section we may easily find the order b correction to the linearized 
solution of the equations of motion. For the metric perturbation hij it takes the form: 

h{i} uu = -2&a (0)n e 4r + 2brq u q v a {0)u e 2r + -brqlq 2 a u{0) + h {lK2)uu e~ 2r + ... 

h{i}uv = ^{i}(2)™e~ 2r ( 52 ) 
h{i} m = h{i}(2) vv e~ 2r 

where we added the curly bracket {1} to indicate that this is the first-order correction. Notice 
that there are no (non-logarithmic) terms at order exp(2r). This corresponds to the fact that we 
set any order b corrections to the source to zero, in agreement with the analysis of [39]. However 
we did allow for a correction to the normalizable terms, which we indicated with a term hny^ij- 
In general these corrections would be unknown but in three-dimensional gravity the equations of 
motion furthermore dictate that: 

h{i}(2)uv = -gbq 2 q%a(fy v + 2ba^^ v 

1vh{i}(2)uv - Quh{i}( 2 )vv = ( 53 ) 

which will become the first-order corrections to the Ward identities of the dual field theory. The 
appearance of the vector field components in /i{i}(2)™ already indicates that the Ward identities 
will change at this order. 

For the vector field the corrections take the form: 

, 1 g 7 2r 2t 

a {i}« = b a q u-qvn(o)w r e + 0{i}(4) u e + . . . 

1 i 2 2 2r (54) 

a {i}v = &fr{o}(2)w _ br(-q u q v h{ y(2)w ~ -^QuQvh(o)w)e r + a{i}(4)»je r + ■■■ 

The term 6exp(2r) in A u , from which all the corrections we wrote above originate, is in our 
conventions a change in the background vector field and therefore does not appear in the above 
expansion. 

5.2 Holographic renormalization 

The computation of the counterterm action is a straightforward implementation of the general 
methods discussed in [39\ 140] . Using Mathematica we found that the counterterm action to first 
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order in b is given by: 

S { i }ct = \j (fx^fb 1 (h ij & j d k a k - rh ij U 1 &d k a k ) . (55) 

Here we replaced the parameter b with a one-form field b{dx l = bexp(2r)du which allowed us to 
write the counterterm action in a covariant form. As usual, indices are raised and lowered and 
covariant derivatives are taken using the background metric 7jj. 

5.2.1 Scheme dependence 

The counterterm action (|55p is not unique, since we may freely add certain scheme-dependent finite 
counterterms as well. For example, a term of the form: 

J d 2 x^^b i h ij n~ / d j d k a k (56) 

is finite as r — > oo and can therefore in principle be freely added to (|55p . Although doing so will 
modify the correlation functions at this order in b by nothing more than contact terms, at higher 
orders in b the addition of such terms will affect the form of the counterterm action and also the 
Ward identities. In order to make our results meaningful at every order in b we therefore have to 
properly identify the scheme dependence. In the next few paragraphs we explain our approach. 

First of all, at every order certain finite counterterms are necessary to maintain general co- 
variance. Since the underlying problem is generally covariant, our counterterms are in principle 
the expansion of covariant counterterms. For example, a volume counterterm leads to a specific 
expansion: 

Tdet^ _> ^det^l + \^h l3 + i( 7 lj M 2 - |7%fc7 fe %i + • • • ) (57) 

as — > -jij + hij. Even though every single term on the right-hand side appears covariant, it is 
only the specific linear combination given above that sums to a completely covariant formal To find 
the relations between the counterterms such that general covariance is retained, we need to verify 
that we are not accidentally breaking diffeomorphism invariance. This can be easily verified by 
demanding that the energy-momentum tensor of the boundary theory is exactly conserved. As we 
will explain in the next subsection, the correct energy-momentum tensor for the deformed boundary 
theory is the tensor Jy defined in equation (foTj) . Demanding then that 

PJij = (58) 
4 Such relations between counterterms are not exceptional and also occur in e.g. the non-linear O(N) model. 
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to any order in b fixes a large amount of scheme dependence in the counterterms. Notice that Jij 
already takes into account the transformation properties of the null vector b under diffeomorphisms 
- indeed, diffeomorphisms in the deformed theory can of course only be unbroken if we allow 
ourselves to transform b as well. 

After having imposed general covariance we are left with a finite number of finite counterterms 
at every order in b. Together with perturbative field redefinitions these reflect the true scheme 
dependence that arises in any covariant renormalization method. In our computations we have 
identified all those counterterms and we added them with arbitrary coefficients. This allowed us to 
fully understand the scheme dependence at every order in b. We will indicate any relevant scheme 
dependence in our results in what follows. 

5.3 One-point functions 

From the first variation of the renormalized on-shell action we obtain the following one-point 
functions: 

1 SS 

(O u ) = = -4(a {0}(4)u + a{i}(4) u ) + (contact terms) 

V ~9(o) 0a (o)v 

1 SS 

(O v ) = = -4(a {0}(4 ) 1 , + a { i }{4 )„) + (c.t.) 

\/-9(o) da {o)u 

(T uu ) = ^1 — -— - = 47r(/t {0} ( 2 ) uu + &/i{i}(2W u - 26a {0 }(4)« + ( c -t-)) (59) 



y-#(o) 




47T 


SS 






4-ir 


6S 



( Tuv ^ = i—^snh, ; = ~ 47r ( h i°}( 2 ) uv + h {m>v) 

y — 9(o) on (o)uv 

(T vv ) = — = 47T (/l{0}(2)ra + hli\(2) vv ) ■ 

\/—9(0) on (0)uu 

Let us briefly explain the structure of these results. First of all, to avoid clutter we have chosen 
not to write the various contact terms since they are largely scheme-dependent. (We will reinstate 
the relevant contact terms in the Ward identities below.) Apart from the contact terms, then, the 
results at order b° are precisely those of equation (|20p written out in components. At order b we of 
course see that the coefficients of the normalizable modes /i{o}(2)ij an d (J{o}(2)i get corrected with 
the hfnmij an d the a{i}(4)j, as expected from the general arguments of [39]. We also find a new 
non-contact term appearing in the one-point function of T uu . This term is indicitative of some new 
structure in the two-point functions, in particular the two-point function {T uu O u ) will be non-zero 
at this order. Apart from the Ward identities discussed below, this non-contact term is the main 
nontrivial result of a careful holographic renormalization procedure at this order. We are not aware 
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of any other (possibly faster) way of obtaining it. 



5.4 Ward identities and Callan-Symanzik equation 

Equations (|53p plus their order 6° counterparts (|15p can now be used to find the Ward identities 
in the deformed theory. We first of all obtain that: 

= q u {T vv ) + q v {T uv ) 

(60) 

= q v (T uu - 2nbO u ) + q u {T uv - 2irbO v ) . 

In agreement with general expectations discussed in [5] and in the appendix to this paper, we 
see that the energy-momentum tensor is no longer conserved. Instead, the correct energy- 
momentum tensor is the non-symmetric combination 

Jij = Tij - IvO.bj (61) 

which satisfies d l Jy = and therefore gives the expected conserved currents associated to transla- 
tion invariance. 

The deformation parameter b also breaks the scale invariance of the dual theory. Indeed, the 
classical beta function for Oi leads to a modified trace Ward identity for the energy-momentum 
tensor, 

{if) ~ 4irb(O v ) =A = 2tt j R[/i (0) ] - cibq 2 a qla {0)v + c 2 6g u ^a (0)n , (62) 

with scheme-dependent coefficients c± and C2- In the appendix we discuss how the left-hand side 
is indeed the expected trace Ward identity for a deformation with a vector operator with scaling 
dimension A = 3 in d = 2. 

The Callan-Symanzik equation takes the form: 

d 2 x( 2 (g<o)ij + ho)ij)jr- — h 2 ( b i + a (o)i)r-— ) s ™ + B = ( 63 ) 

v on (0)ij oa (0)i / 

with the integrated conformal anomaly 

^ = 2 J ^ 2 9(( a (0)j;9M^ a (0)t; + a (o)uQuqf,a( )u) ~ bh^ vv qlq v (q u a^) v - q v a.( )u) + ^[^(0)]^ a (0)u) • 

(64) 

The new terms in the conformal anomaly are scheme- independent. From equation (|63p we read 
off that the sources still scale with the same scaling dimension as for b = so there are no beta 
functions at this order. 

We will present a more extended discussion of the Ward identities and the Callan-Symanzik 
equation for finite b in the next section. 
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5.5 Two-point functions 

As before, the two-point functions can be found by expressing the normalizable modes and 
h(2)ij m terms of the sources tt( )i an d fyo)jj> substituting in (|5T?|) and taking another derivative 
with respect to the sources. Using the fact that we can solve the linearized equations of motion 
completely for finite 6 (see section 0] for the asymptotic expansion and section below for the 
regularity condition) we find the following first-order expressions for the normalizable modes: 

«{o}(4)« + a{i}(4)« = g?«9«(ln(?„g„) - ln(2) - 1 + 27) (o( )„ - ^bh^ vv ) 



+ 



(65) 



1 1 ^ 

a{o}(4)^ + a{i}(4)v = g?«9«(M9w9«) ~ ln ( 2 ) ~ 1 + 2 7) («(o)« - &tyo)u« + h {0)uu ) 
o z q u 

+ -^ > 1u ( ll a (0)v ~ -^bqlqlh {0 ) vv . 

The hriy(2)ij are again completely determined in terms of the sources by the constraints ([53 
Altogether this leads to the following new non-vanishing two-point functions: 

d 2 q , , , . s .■„.„ 36 



u^v 2 



(T uu O u ) = -ibir J -^^qlq v log(q u q v 

(T vv O v ) = ibit [ -pLqt \og{q u q v - ie)e^ x = ^ (66) 
J (2vr) z uv 

f d 2 q o , „„ „ 36 

{T uv O v } = -ibir J -j—^q u qtAog{quqv 



I 9 4' 

whereas the other two-point functions are unchanged and we again ignored various contact terms. 
In terms of the currents Jij we find the following complete set of non-zero two-point functions: 

(JuuJuu) — 7~ {&u@u} — ~ A 9 

U ITTU^V* 

(JvvJw) = (O v O v ) = 1 T-T (67) 

U Z V^ uv 

The only new correlation functions are those on the last line. The first of these indicates that, 
although (J uv J U v) vanishes up to a contact term, the descendant operator d v J uv is not a trivial 
operator. In the analysis of [IB] it was claimed that in a two-dimensional theory with non-relativistic 
(or 'chiral') scale invariance, at least under certain assumptions, d v J uv had to be a trivial operator 
and the conservation equation d % J% v would then dictate that d u J vv vanished as well. This would 
lead to the existence of an antiholomorphic current J vv which together with J uu would lead to two 
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infinite-dimensional symmetry algebras. The last line of (|67p however explicitly indicates that this 
scenario is not realized in the deformed theory already to first order in 6£ 



6 Analysis for finite b 

The next step in the analysis is to increase the order in b in which we are working. Using Mathe- 
matica we have worked out the full holographic analsyis up to and including order 6 4 . Our reason 
for working to this order is that all possible contractions of hi with at most two out of the fields 
hij and a, are exhausted at order b 4 , with the final structure given by: 

J (Fx^VVb^hijhkt (68) 

and what remains at higher orders are infinite series expansions in the scalar operator b 2 q^. (The 
above counterterm is indeed necessary at order 6 4 , its coefficient happens to be —1/12.) 

The analysis at order 6 4 leads to several interesting observations. For example, up to order 6 4 
we find from equation ([76]) below that there is no explicit term of the form b u a^ v in the one-point 
function of T uv , even though such a term would be allowed by conformal and Lorentz symmetry. 
Although our analysis does not exclude the appearance of such a term at order 6 5 , for example with 
an extra factor 6 4 g 4 , its absence at the lowest two orders leads us to believe that this term will be 
absent to all orders. Using this and similar observations, we will in this section not only present 
the results of the perturbative analysis but also conjecture a consistent extrapolation of our results 
that should be valid to all orders in b. The equations in this section represent the results of this 
extrapolation, whilst the truncation of the results to order 6 4 is based on explicit computations. 

The discussion below in principle follows the pattern of sections [3] and [5j We have however 
already extensively discussed the asymptotic solution to the linearized equations of motion in section 
H] and we will therefore proceed with a discussion of the holographic renormalization procedure. 

In the previous section we indicated the order b correction to for example h^ij by /i{i}(2)ij and 
correspondingly introduced /i{o}(2)«j to indicate the order b° result. In this section we will drop the 
curly brackets and simply write h( 2 )ij to denote the all-order result, so: 

h(2)ij = h{o}(2)ij + h {l}(2)ij + h{2}(2)ij + • • • (69) 

and we employ a similar change of notation for other coefficients. 

For finite b we will see below that the theory becomes non-local and the spectrum of the non-relativistic dilation 
operator becomes continuous. This explicitly violates the assumptions of [16] so the results of that paper will then 
automatically no longer apply. See also [62] for a related discussion. 
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6.1 Holographic renormalization 

As we mentioned above, the procedure for obtaining the counterterm action in the presence of 
irrelevant deformations was explained in our earlier work [391 HO] . With the aid of Mathematica 
we obtained the full counterterm action to order 5 4 , using the methods explained in subsection 15.21 
to ensure covariance of the counterterm action. Furthermore, at each order we have determined all 
the possible finite local counterterms in order to be able to parametrize all the scheme dependence 
in our results. 

The full form of the counterterm action is rather complicated and not too illuminating so we will 
not present it here. Let us however mention that, starting at order 6 2 , we also need double-trace 
counterterms. The full double-trace part of the counterterm action to order 6 4 takes the form: 




where the Cj are scheme-dependent coefficients and Ojj and irt are the renormalized momenta con- 
jugate to hij and cij, respectively. These are defined as [30] : 

6 S 

0* 3 = TT— (Share + Sct.st) VT* = — (Share + Sct.st) (71) 

ohij dai 
with 5 c t, s t the single-trace piece of the counterterm action. Notice that most of these counterterms 
are logarithmic in nature, i.e. they involve explicit powers of r, leading to operator mixing between 
single- and double-trace operators in the Callan-Symanzik equation below. There are no higher- 
trace counterterms since the background conjugate momenta vanish identically and we only expand 
the action up to second order in the fluctuations. 

As explained in [3U] , logarithmic double-trace counterterms generally force one to perform field 
redefinitions in order to retain a correct variational principle. Indeed, for example at order b 2 we 
find terms in the variation of the on-shell action of the form: 

<5Sren = / d 2 X[ . . . + h(2) uu 8h(Q) vv + hm m 8h(Q\ uu - 2h(2) uv 5h(Q) uv + ... 

(72) 

... - 6 2 (1 - 2c )/l( 2 ) w 5/i(2) OT +...)■ 

From the order b° metric variations on the first line we find that h^)w 1S a conjugate momentum so 
its variation should not appear in an action with an appropriately defined variational principle. In 
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order to repair the variational principle at order b 2 we may either choose the scheme where cq = 1/2 
or alternatively we redefine the sources to be: 

1. 



h(o)uv -> h(o)uv - 2 &2 ^ ~ 2c o) h (2)w (73) 

so the variation 8h(p\ m cancels. The choice between the different prescriptions amounts to adding 
a finite double-trace counterterm to the theory, see [lOj for details. 

In this paper we will opt for the simplest solution and set Co = 1/2 so no redefinition of the 
source terms is necessary. At higher orders we can continue this procedure and fix the finite terms 
in the multi-trace counterterm action precisely such that no field redefinitions are needed. This 
leads to the simple prescription that the sources are always the terms of order exp(2r) in both Oj 
and in hij. 

6.2 One-point functions 

Let us now present the one-point functions, defined as the first variation of the renormalized on-shell 
action. We find the following one-point functions for the vector operator: 

W = -^^7^ = -4[1 + ~ ^7-b\t + . ..]a {4)u 

5 2 

+ -^> 2( it a v [l + -K^ll + ■ ■ ■] a>(4)v + (contact terms) (74) 

<a> = —±=J!L = -4[i + h\i - hy v + . ..]a (4)v + (c.t.) . 

^-9(0) °a(0)u 4 8 

The contact terms here are scheme-dependent and we have omitted them. Although we have 
computed the on-shell action only to order b 4 , we expect that at higher orders in b no completely 
new terms will appear in the one-point functions. On the other hand, it seems likely that the 
expansions in b 2 q 2 that we enclosed in square brackets in equation (|74p will continue indefinitely 
to form a non-trivial function of b 2 q 2 . We notice that this function can in fact be removed in 
perturbation theory by a redefinition of the sources 

O(0)i -> [1 + #b 2 q 2 + #bY v + ■■■} o (0 )i (75) 

for some appropriate numerical coefficients we denoted with a # here. Such redefinitions are local 
to any finite order in b and are therefore an allowed part of the usual scheme dependence in quantum 
field theory, complementary to the addition of finite local counterterms in order to shift the contact 
terms in the correlators. 
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For the one-point functions of the energy- momentum tensor we find that: 

(T uu ) = ~^= x f S = 4vr(/i {2 ) mt - 26 [l + ]rb 2 q 2 + . . . ] a (4)u - ^b 3 [l + . . . ] q u q v a {A)v + (at/ 
y/9(0) on (0)vv \ y 

(T uv ) = -^t^~ = -47T ( [1 + -b\l - \b\i + . ..}h (2)uv + (at.)) 

(T™) = ^T7^ = 47T f [1+ ^ - -W 4 + • • .]/»(2) W + (Ct: 

(76) 

where we again omitted the contact terms and wrote in square brackets the possibly infinite ex- 
pansions in b 2 q 2 . Notice that we did not find any b 2 q 2 terms multiplying h^)uu m {T U u) U P to and 
including order b 4 . Furthermore, the expansions in b 2 q 2 multiplying h( 2 ) U v an d h^)w m {T U v) an d 
(T vv ) are identical at least to order b 4 q 4 . We will again assume that this is no accident and that 
this effect persists to all orders, even though we cannot currently offer a detailed explanation. 

6.3 Ward identities 

The equations (|4ip combined with the redefinitions (|48p ~ (|5ip result in constraints for the h^)ij- 
Using then (|76p and (|74p we can rewrite these constraints in terms of the one-point functions of 
the dual operators. This leads to: 

= q u {T vv ) + q v {T uv ) 

(77) 

= q v {T uu - 2irbO u ) + q u {T uv - 2nbO v ) . 



The correct energy-momentum tensor is therefore again which we recall was defined in (|6ip to 
be equal to Jij = Tij — 2irOibj. This tensor is the conserved energy-momentum tensor, 

9 l J 4J =0, (78) 

which we now claim will hold to all orders in b. (As discussed in section [5j this conservation 
equation is essentially a consequence of our covariant renormalization scheme.) 

The trace Ward identity again accurately reflects the breaking of scale invariance. With respect 
to the analysis at order b 1 we find that it is augmented only by contact terms: 



T\ - 4irbO v = A = 2n[l + #b 2 q 2 + #b\ A v + . ..]%„)] 

(79) 



+ 6[# + #b 2 q 2 v + . . .]q 2 u q 2 a {0)v 



+ b[# + #b 2 q 2 + ...] 
+ b 2 [# + #b 2 q 2 + ...]q 2 u q 2 h {0)vv 
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where various (scheme-dependent) numerical coefficients are denoted by As we mentioned 
before, in the appendix to this paper we explain why the left-hand side is the expected trace Ward 
identity for a deformation with a vector operator with scaling dimension A = 3 in d = 2. The 
right-hand side is the unintegrated conformal anomaly to first order in the sources. 

As explained in detail in the appendix, we can combine the above equations to find that the 
current Ri, defined as: 

Ridx 1 = J uu udu + J vu udv (80) 

is conserved up to contact terms, 

0*2^ + ^/4 = 0, (81) 

reflecting the non-relativistic scale invariance of the theory for finite b. In fact, we may rewrite (|79p 
as: 

Jvu = \A (82) 

which demonstrates that J vu is a trivial operator, i.e. all its correlation functions reduce to contact 
terms. Since d l Jij = 0, we find that the operator J uu satisfies d v J uu = and therefore generates 
an infinite-dimensional Virasoro symmetry for finite b as well. This 'left-moving' Virasoro algebra 
is thus unaffected by the deformation. 

6.4 Callan-Symanzik equation 

The methods for obtaining the Callan-Symanzik equation in the presence of double-trace counter- 
terms were extensively discussed in [JO]. We find that it takes the form: 



J d 2 x (2fl>( 0)ij + 2h {0)ij 

+ J d 2 x (2b + 2[i + ...]a {0)u + h 3 [i + ...]d u d v h {0)vv + #b 4 [i + ...]d u d, 

+ J d 2 x (2[l + ...]a {0)v + ^b 3 [l + ...]d 2 v 
+ J d 2 x(-\b 2 [l + ...]R[h {Q) ]+b 2 [l + ...]dlh 



3„ 
a (0)v 



5S r , 



a 2. \ fl£rcn 



(83) 



(0)vv 



+ #6 3 [1 + . . . ] d u d 3 a {0)u + #6 3 [1 + . . . ] d 2 u d 2 a {0)u ) 



J 5h(py 



+ f d 2 x2b 3 5 x Sien x f ren =-B 
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Here the dots represent infinite series in b 2 q 2 and the ft represent scheme-dependent numbers. The 
conformal anomaly B can be conveniently represented as: 



B = s l Bs 



with s the vector of local sources, 



and the matrix B given by: 



B 



-' = 1 


°(0)u> a {0)v 


h(o) vv ,R[h( )\) 


2 u u u v 


ftb 2 dldt 


-ibd 2 d 2 ifea 2 \ 

2 uu u u v 2 uu v 





2 u u u v 


-\bdld v #b 3 d u d* 








\b 2 dld v ftb^l 


V o 





0{b e ) J 



(84) 



(85) 



(86) 



where, although we have not written it explicitly here, each term is expected to receive an infinite 
series of corrections in b 2 q 2 . Notice that all the terms in B have the correct weight and tensor 
structure so that B is a scalar of dimension 2. 

Let us now explain these equations a bit further. First of all, we recall from [40] that the final 
term in (|83|) . which is quadratic in functional derivatives, represents operator mixing between single- 
and double-trace operators at large N. From this term we conclude that the deforming operator 
O v mixes with the double-trace operator O v T vv at order 6 3 . Furthermore, the new coefficients 
multiplying the other functional derivatives of the on-shell action in ()83f) represent the anomalous 
scaling dimensions of the corresponding operators. We see that all the anomalous scaling dimensions 
in (|83p involve derivatives and therefore vanish for constant sources. This is in agreement with the 
all-order analysis of [HE], although in those papers the (multi-trace) operator mixing was not taken 
into account. 

Of course, the conformal anomaly and the beta functions involve many terms and there is no 
clear structural pattern. (The one exception is that h^ uv and h(Q\ uu only appear in the linearized 
curvature R[h^], which is presumably a consequence of our covariant renormalization scheme.) 
Such a plethora of terms is actually to be expected, since an irrelevant deformation generally leads to 
the addition of all possible counterterms consistent with the residual symmetries. Whenever these 
counterterms are associated with logarithmic divergences, a corresponding term in the Callan- 
Symanzik equation appears and this leads to the multitude of terms we observe in the above 
equations. 
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At higher orders in b we expect no new tensor structures to appear in the Callan-Symanzik 
equation. Instead, we expect that all the terms in the equations (|83p and (|84p will get corrected 
with infinite series in the scalar operator b 2 q%. For finite b, such series will generally resum to 
terms that are non-local in the v direction. This non-locality in the Callan-Symanzik equation is 
an unavoidable consequence of the finite irrelevant deformation. It signals that for finite values of 
b the w-dependence of our results will be scheme-dependent and essentially devoid of any physical 
information. Fortunately, as explained already in [5], our deformation retains locality in the u 
direction. The u-dependent part of our results is therefore scheme-independent and representative 
of the actual dynamics of the deformed field theory. 

6.5 Two-point functions 

The two-point functions can be computed by expressing the normalizable modes and h^)ij 
in terms of the sources O(o)j and /i(o)ij by demanding regularity of the solution in the interior of 
the spacetime. In our case regularity in the interior can be easily imposed by going back to the 
solution written in terms of Bessel functions (|34p and expanding it around r = — oo. In terms of the 
integration constants c±\ and c±g introduced in (|34p we find that the necessary conditions become 
c+i = c_i and c+g = c_g. In terms of the ^(0)1 an d ^(4)i integration constants we then find that: 

= ^2^9^(0)!) ip(4)v = 7^2X1 V>(o)« ( 87 ) 



where we introduced: 



2 2 -^ 2 (3 + S9 )r(i- S9 ) . 
X9 = j^p^j g tt (gwb - «) 

2-*i6 2 r(i - Sl ) 2 . ! 
r(2 + Si) 

to simplify the expressions below. These expressions, combined with the expressions (|4ip which 
express the X(2)ij m terms of the X(o)iji can be substituted in our definition for the sources (|48p - (|5ip 
to express all the normalizable modes in terms of the sources. We can then take a further functional 
derivative of the one-point functions (I74p and (I76p to arrive at the two-point functions given below. 

For the two-point functions of the vector field we find it to be advantageous to introduce a new 
operator O u which will replace O u , defined as: 

6u = [1 - ^f q l - ^bV v + ■ ■ ■ ] (0 U + \b\ u q v [1 + b\l + . . . ] a) (89) 
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which re-diagonalizes the two-point functions in the (O u , O v ) sector and also introduces a convenient 
normalization for O u . Furthermore, we also define a new normalization for the other operators by 
introducing: 



where the factor on the right-hand side is the square root of the factor [1 + 6 2 g^/4 — b 4 q 2 /8 + . . .] 
appearing in various one-point functions in flT4f) and (1761) , 

With these conventions we find the following correlation functions for two vector operators: 



where on the right-hand side of these equations we omitted scheme-dependent contact terms. Just 
as for the Callan-Symanzik equation at finite 6, these contact terms will generally involve infinite 
series in b 2 q 2 , so for finite b they become arbitrary functions of v 2 /b 2 . (Furthermore, analyticity 
at the origin of momentum space dictates that the corresponding position-space expressions are 
integrable in the v direction.) This is again precisely the infinite scheme dependence and non-locality 
typically associated with a finite irrelevant deformation. On the other hand, all of these infinitely 
many contact terms are at most polynomial in q u and therefore the position-space expressions will 
indeed be contact terms in the u direction. 

Two-point functions involving one (new) energy-momentum and one vector operators take the 
form: 



x = x[i-hy v + ^-b\t + ...] 




(90) 




(91) 



(O u O v ) 











(92) 
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Finally, the two-point functions involving only the energy-momentum tensor become: 

d 2 Q in-x <?« _ 12vr 



(27r) 2 (q u q v - ie) U 



1 



(JuuJuu) — 47T i 

{JuuJuv} — 
(JuuJvv) — 

If J \ - _0^j [ d2(1 P fr* QUQV (93) 

{J UV J UV ) - M i J {l _ Xi) 
d 2 q " 2 



{JvvJvv) — 27T i 



2 ,' / _ 1 J** 



(2vr) 2 (l-xi) 

2- /" ^ 2{ ? Jqx it 



(2vr) 2 (9u5«-*e)(l -xi) 
where we again suppressed contact terms on the right-hand sides. One may explicitly verify that 
the Ward identities discussed above are satisfied. 

Performing the inverse Fourier transforms in the above correlation functions seems quite com- 
plicated, although it is not hard to see that for small q u and fixed q v , so at large separation in the 
u direction, all correlation functions reduce to their AdS3 counterparts. We can also try to scale 
out the u-dependence to find a mixed position-momentum space expressions. In this way we find 
that for example the first equation in (]9ip becomes: 



(6 u (u, q v )6 u {u\ -q v )) ~ \u - uT 3 '- 1 ~ h2( 9 n4 (94) 

The last factor can be removed by shifting O u — > O u + ^VSJuu/i^b) although in that case the 
b — > limit of the correlator is no longer smooth. From a non-relativistic viewpoint q v becomes a 
quantum number (the particle number in condensed matter applications [H[2]) and the powers of 
u are the non-relativistic dilatation weights. We see that these weights have become dependent on 
q v , so for an uncompactified v direction the spectrum of the non-relativistic dilatation operator is 
continuous. For applications to condensed matter physics one should demand a discrete spectrum. 
One method for obtaining this spectrum is to compactify the lightlike v direction. A discussion of 
the subtleties associated with this DLCQ can be found in [36]. 



7 Conclusions 

We have holographically computed the two-point functions of the deforming vector operator and the 
energy-momentum tensor for the field theory dual to the Schrodinger backgrounds using a consistent 
irrelevant deformation picture. We discussed in detail the non-relativistic Ward identities satisfied 
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by the correlation functions. The holomorphicity of the operator J uu in particular indicates that the 
deformation keeps the 'left-moving' symmetries of the original CFT intact, which is in agreement 
with the viewpoint that the deformation is exactly marginal from a non-relativistic point of view. 
The right-moving Virasoro is however broken already at first order in b. It is also interesting to 
note that the Callan-Symanzik equation demonstrates mixing between single- and double-trace 
operators. This leads one to conclude that the question whether the deformation is of single- or 
multi-trace form becomes renormalization scheme dependent. It would be interesting to see if our 
results could be matched to computations in dipole theories. 

One very valuable general lesson arising from our analysis is that, to any finite order in any 
irrelevant deformation parameter, the sources are always the term of order exp(— A_r) in the 
radial expansion of the fields, with A_ being the appropriate scaling dimension for the sources in 
the absence of the deformation. In other words, the location of the sources in the radial expansion of 
the fields is unmodified to any finite order in perturbation theory. Only when we consider all-order 
effects in the deformation parameter do we observe a modification of the location of the source 
terms. Interestingly, in such cases the analysis becomes very subtle as the sources generally do not 
correspond to the leading terms of the radial expansion of the bulk fields. We demonstrated the 
resolution of these issues explicitly in section HI 

We would like to stress that the proper identification of the source terms lies at the heart of the 
problem of defining a proper holographic dictionary for these non-AlAdS spacetimes. Notice that 
this is in essence a gravity question, since the correct identification of the sources also determines the 
proper boundary conditions for the bulk dynamics. This is an essential step in order to for example 
define conserved charges, the notion of 'asymptotically Schrodinger spacetimes' and the Schrodinger 
equivalent of 'normalizable modes'. In this light it is interesting to note that, if one expresses these 
proper boundary conditions in terms of the induced bulk fields at finite r, the presence of multi- 
trace counterterms dictates that these boundary conditions will involve the conjugate momenta of 
the fields as well. Dictating simple failoff conditions for the metric, which happens frequently in for 
example an asymptotic symmetry group analysis, is therefore not sufficiently general for spacetimes 
with complicated asymptotics. 

Our results come with a significant scheme dependence which indicates a certain dependence 
on the UV completion of the boundary theory even at the level of classical bulk gravity. It would 
be interesting to see whether there perhaps exists a 'holographic renormalization scheme' which 
would be natural from the perspective of the gravity dual. In this light it is also important to 
realize that at large N one can perform additional redefinitions of the sources which correspond 
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to e.g. multi-trace deformations of the theory. We have fixed this freedom by demanding that the 
sources are the terms at order exp(— A_r) in the radial expansion of the fields, at least to every 
finite order in perturbation theory. It may however be worthwhile to further explore this freedom. 

Our analysis paves the way for the computation of re-point correlation functions in three di- 
mensions. Furthermore, there should be no conceptual difficulties in the extraction of correlation 
functions from higher-dimensional Schrodinger spacetimes or from black holes with Schrodinger 
asymptotics [341 [35| [36] . It would also be very interesting to generalize our results to other space- 
times which are not asymptotically AdS, for example the Lifshitz spacetimes or the spacetimes of 
[63"1 [64] dual to non-commutative gauge theories. 
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A Vector deformations in quantum field theory 

Consider the partition function Z of a (i-dimensional conformal field theory in a background metric 
gij deformed by a source bi for a vector operator Oi with dimension A. In path-integral notation 
the partition function can be abstractly written as: 

z [9ij, bi] = J D§ exp ^ - iS[gij, *] - i J d d x^fgg %i biOj + (ie insertions)^ . (95) 

In this appendix we discuss the modification of the field theory Ward identities that this deformation 
brings about. The ie insertions are important for obtaining the correct correlation functions in 
Lorentzian signature but we will ignore them here. 

In the conventions imposed by the holographic analysis the energy-momentum tensor is defined 
to be the operator dual to the boundary metric, more precisely 



T ij = - 4 "_ r Z\g ih bA = T ij - 27r(g^b k O k - b i O j - WO 1 ) (96) 

^gSg i: 
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with T %3 the energy-momentum tensor before the deformation. Notice that we hold the one-form 
hi fixed while varying the metric. 

Under a diffeomorphism along a vector field Q both the metric gij and the vector field bi 
transform by their Lie derivatives: 

Sgij = V,Q + VjC* SIh = CVjbi + bjViC (97) 

and in the absence of a diffeomorphism anomaly in the path integral measure we should find that 
the partition function is invariant, 

5c_Z = , (98) 

leading to 

V*(2ij - 2TrOibj) + 27r(V i i» i )O i = . (99) 

We will be working in flat space and with a constant 6, so Vjfoj = 0. In that case we find that the 
conserved current associated to translation invariance is given by: 

J i:j = T i:j - 2nOibj d^ij = (100) 

This is indeed confirmed by the holographic analysis in the main text. 

Since the new energy-momentum tensor Jij is not symmetric, Lorentz invariance is manifestly 
broken, even if dibj = 0. For example in d = 2 a boost can be implemented by a diffeomorphism 
with £ 3 ' = e\x k . In that case we find that 

d\T i3 t\x k - 2irb 3 e j k x k Oi) + 2irb J e 3 k O k = , (101) 

which does not integrate to a conservation equation. 

Let us now discuss Weyl covariance to first order in b. A local scale transformation is imple- 
mented by sending g^j —> 2Xgij and 6j — > A(A — d+\)bi with A the infinitesimal scale transformation 
parameter. Assuming that the transformation of the path integral measure results in a local con- 
formal anomaly A, the associated Ward identity would become: 

T\ - 2vr(A - d + 1)6*0* + A[gij] = . (102) 

This equation is valid to first order in b. At higher orders one generally finds that the classical 
scaling dimensions get corrections, so new one-point functions may appear and the classical scaling 
dimension A becomes a matrix of anomalous dimensions. Similarly, the conformal anomaly A may 
receive 6-dependent corrections. 
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We may combine (|102j) with the diffeomorphism Ward identity (jlOOp to find that: 

d i (T ij x j - 2-KXjVOi) - 2^(A - d)b t O l + A[g tj ] = . (103) 

Now let us specialize to the case where d = 2 and b{ is null, bidx 1 = bdu. For such a vector one 
verifies that bje J k = — bt (in conventions where e uv = +1) and therefore we can combine (I103|) and 
([TOT]) to find that: 

d l R t + l -A[g l3 ] = 0, (104) 

with 

ifc = Tij(x j - (A - d)e\ x fc ) - 2irbjX j (A -d+ 1)0; , (105) 

which is the current associated to the non-relativistic scale invariance of the theory. In the main 
text we have A = 3 and d = 2. In that case the current Ri becomes: 

Ru = JuuU Ry = J vu u (106) 

In fact, precisely for these values for A and d we see from equation (|102|) that 

Jvu = -^by] ( 107 ) 

and J vu is therefore 'trivial' in the sense that its correlation functions reduce to contact terms. The 
conservation equation d 1 3{ u therefore dictates that, 

d v J U u = (108) 

up to contact terms. We thus immediately see that we can build infinitely many conserved currents 
of the form: 

Juuu k (109) 
for some integer k. In other words, the left-moving Virasoro algebra is left intact by the deformation. 
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